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Spin-state transitions, observed in many transition metal compounds containing Co''"'" and Fe'^"'", 
may occur with the change of temperature, pressure, but also with doping, in which case the 
competition of single-site effects and kinetic energy of doped carriers can favor a change in the spin 
state. We consider this situation in a simple model, formally resembling that used for manganites 
m Ref. H. Based on such a model, we predict the possibility of a jump-like change in the number 
of Co^+ ions undergoing spin-state transition caused by hole doping. A tendency to the electronic 
phase separation within a wide doping range is demonstrated. Phase diagrams with the regions of 
phase separation are constructed at different values of the characteristic parameters of the model. 
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I. INTRODUCTION 

Interplay of different degrees of freedom and different 
types of ordering is a very important ingredient in de- 
termining the properties of strongly correlated electron 
systems. Especially interesting these effects become in 
doped systems. Typical for this case is the tendency to 
phase separation and formation of inhomogeneous states. 
It can take different forms: formation of isolated polarons 
or small clusters (modification of a particular ordering by 
doped charge carriers and the trapping of charge carriers 
in a distorted region) or particular textures, e.g. stripes. 
Such phase separation can play a very important role in 
many phenomena, such as colossal magnetoresistance in 
manganites^, and probably also in high-Tc superconduc- 
tors - although their role in the latter is still a matter of 
hot debate. 

The most common and best known case is the dop- 
ing of antiferromagnetic insulators, with the formation 
of ferromagnetic droplets ("ferrous") or charged antifer- 
romagnetic domain walls (stripes)2ii2ii^. We have recently 
shown that, similarly, the interplay of kinetic energy of 
doped holes with the orbital structure can give rise to a 
novel mechanism of phase separationSii^. 

A special interesting group of phenomena is met in 
systems where the respective ions can exist in different 
multiplet states. Typical examples are the compounds 
containing Co'^"'" (or sometimes Fe^"*"), which can exist in 
a low-spin (LS) state with S=0 (^23)1 intermediate-spin 
(IS) state 5 = 1 (^igCp, and high-spin (HS) state (^363) 
with 5 = 2, see e.g. Ref. ItI. Close proximity in energy of 
these states can lead to a special type of transition (or 
crossover): spin-state transition (SST), typical example 
being LaCoOa^'^'^'iS. Also spin-state ordering is possi- 
ble^iiS. Thus, these systems, in addition to quite com- 
mon charge, orbital, and spin degrees of freedom with the 
possibility of respective orderings, have an "extra dimen- 



sion": the possibility of spin-state (or, in other words, 
multiplet) transitions. Correspondingly, if doping of ma- 
terials like manganites can cause phase separation due 
to an interplay of the motion (kinetic energy) of doped 
holes with the underlying magnetic and orbital structure, 
in systems with SST like cobaltites one can expect simi- 
lar phenomena due to an interplay with the spin state of 
the matrix. The common mechanisms causing the phase 
separation manifest themselves in the situation when the 
particular ordering existing in the system hinders the mo- 
tion of doped holes. In these cases, it may be favorable to 
locally modify the type of ordering, facilitating the mo- 
tion of the hole in such distorted region. Thus, holes can 
hardly move on an antiferromagnetic background, which 
was noticed already long ago both for the two-band (dou- 
ble exchange) mode l^^'^"^ and for the single-band (Hub- 
bard) modeli^ii^. At the same time a hole moves freely 
on the ferromagnetic background. As a result, ferromag- 
netic polarons (ferrous) may be formed close to the hole, 
and the gain in kinetic energy of the latter moving on the 
ferromagnetic background exceeds the loss of the mag- 
netic energy 
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Similarly, certain types of orbital ordering suppress 
hole motion, and it may be favorable to modify orbital 
pattern close to a hole, forming orbital polaron and facili- 
tating motion of a hole within it^^d^n^. For systems with 
SST such a role can be played for example by the phe- 
nomenon of a spin blockade^!: if one dopes the material 
with the Co'^"'" in a low-spin state {S = 0) by electrons, 
the ionic state created could be Co^+ in a high-spin state 
(5 = 3/2). In this case, it is evident that it is not pos- 
sible to interchange the states Co'^"'" LS and Co^+ HS 
by moving only one electron: one would end up in the 
"wrong" states Co^+ and Co^+ both in IS states, not in 
the original states (the hopping of an electron can change 
the spin of corresponding states only by ±1/2, whereas 
the spins of the original states differ by 3/2). As a con- 
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sequence, an extra electron can only move in a crystal 
leaving the trace of wrong spin states, which will lead to 
a confinement and localization of this electron, similar to 
the case of the usual Hubbard modeii^. 

One can "repair" this by modifying the spin state 
in the vicinity of a charge carrier (electron or hole), 
and this will finally again lead to a creation of inho- 
mogeneous states and to phase separation. This phe- 
nomenon was actually observed in some cobaltites, e.g. 
in Lai_a;Sra;Co03. There are already many indications 
of phase separation and formation of inhomogcncous 
states in this systen*2^i2^i2ii^i2£, but probably the most 
clear evidence comes from the study of very low doped 
LaCoOa. Magnetic measurements^ have shown that at 
very low doping (< 1% of Sr) the moment per doped hole 
(per Sr) is much bigger than that of only a LS 00**+ with 
S — instead there exist magnetic impurities with 

unusually large spin S = b — 10, which signals that each 
hole is "dressed" by the magnetic cloud due to the pro- 
motion of some of neighboring Co^+ ions to a magnetic 
state. The magnetic scattering, ESR and NMR study 
of such system^^ allowed even determining the size and 
shape of such magnetic clusters formed around doped 
holes. 

It is possible to use different approaches to describe 
theoretically the phenomenon of phase separation. First 
of all, it is the direct numerical investigatioiiS. Or one 
can assume the formation of spin-state polarons, calcu- 
late their energy and check whether and at which condi- 
tions the formation of such polarons can be energetically 
favorable. But the most direct way, by which one usually 
starts, is first to assume the existence of a homogeneous 
state and to check for its stability against phase separa- 
tion. This was the route taken earlier by us for the double 
exchange model^, for the situation close to a charge or- 
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deringsi, for two-component model of manganite: 
for orbital ordering^. If the homogeneous state turns out 
to be unstable, then at the second step one can investi- 
gate particular types of inhomogeneous states, which can 
be formed. In the present paper, we follow this route for 
the doped systems with SST. 



II. SPIN STATES OF COBALT IONS 

Let us list the possible spin states of Co^+ and Co^+ 
ions in a CoOg octahedron, which is a main building 
block of perovskite-like Co-based compounds (we will 
consider below the hole-doped cobaltites, nominally con- 
taining Co'^"'' and Co'*+). The electron configuration of 
Co'^+ ion is 2>cP . It is well known that in the crystal field 
of cubic symmetry, a d-level with the 5-fold orbital de- 
generacy is split into a doubly degenerate eg level and 
a triply degenerate t2g level. In the octahedral coordi- 
nation, the t2g level lies below the eg level. So, a Co^+ 
ion can have three low-energy spin states: low-spin (LS), 
intermediate-spin (IS), and high-spin (HS) states. 

In the LS state {S = 0), all t2g states are occupied 
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Table I: (Color online) Possible electron configurations of Co 
ions and their energies. 



and the eg level is empty. In the IS state {S = 1), there 
are five electrons at the t2g level and one eg electron. 
In the HS state {S — 2), we have four t2g and two eg 
electrons. The corresponding energies of these states are 
Els'' = ^0, E^f^ = Eo + A - Jh, and E^j^p = Eo + 
2 A — 4:Jh, where A is the energy splitting between t2g 
and eg levels and Jh is the Hund's rule coupling constant. 
For the Co^"*" ion (Sd^), there are three similar low- lying 
spin states, corresponding to different distributions of five 
electrons between t2g and eg levels. In the LS state {S — 
1/2), there are five electrons at the t2g level and no 



electrons. In the IS state (S 



3/2), we have four t2g 
In the HS state (5 = 
5/2), the numbers of t2g and eg electrons are equal to 
three and two, respectively. The corresponding energies 



electrons and one Cg electron 



of these states for the Co''"*' ion are E^j^p 



Ei + A-2Jh, and e'-^P = Ei+2A~ GJh- Here, we 
introduced Eq and Ei as some reference energy values for 
Co'^^ and Co'*'*', respectively. As we shall demonstrate 
below, the results do not depend much on the specific 
choice of Eq and Ei. All aforementioned configurations 
of Co ions and their energies are summarized in Tabic [D 
The type of the ground state for a separate Co'^+ or 
Co''^"'" ion depends on the relationship between A and 
Jh. It can be easily seen that at A > 3Jh, the LS is the 
ground state both for Co^+ and Co^""". At 2Jh < A < 
3Jh, Co'^^ still has the LS ground state, whereas for 
Co^''' the HS is more favorable. Eventually, at A < 2Jh, 
the HS state is the lowest in energy for both ions. Hence, 
for isolated cobalt ions, the IS ground state does not 
appear. 

The situation becomes more complicated if there exists 
a charge transfer between cobalt ions. First, note that 
the hopping integrals between the t2g states in cobaltites 
are as a rule much smaller than for the eg states. In 
the treatment below, we ignore the t2g — t2g hopping and 
take into account only the hopping of eg electrons. The 
inclusion of t2g — t2g hoppings will not modify qualita- 
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tive results, introducing only minor numerical changes. 
Second, the states with the number of electrons per ion 
larger than six are unfavorable due to the strong on-site 
Coulomb repulsion. Third, the transitions of electrons 
between the lattice sites corresponding to the changes of 
spin by more than one half are strongly suppressed since 
they involve the simultaneous change of a state for two 
or more electrons. 

As a result, in doped cobaltites there remain only two 
most probable hopping processes: (i) the transitions of 
electrons between the IS Co^+ and LS Co''+ and (ii) tran- 
sitions between the HS Co^+ and IS 00**+. The corre- 
sponding configurations are illustrated in Table [J Thus, 
to facilitate the kinetic energy gain due to the electron 
transfer, one can create a ground state with intermedi- 
ate spins. Such a situation can arise if in the ground 
state for isolated ions, we have either LS Co^+ or HS 
Co'^"'". The former case corresponds to A > 3Jh when 
some of LS Co'^+ can be promoted to the IS state. In the 
latter case corresponding to A < 2Jh, some HS Co^+ 
are promoted to the IS state. In the intermediate situa- 
tion, 2Jh < a < 3Jh, the electron transfer can occur if 
we promote both ions, Co"^"*" and Co'*"'", to some excited 
states. Such double excitations seem to be less probable. 
Below, we first discuss the most realistic case A > 3Jh 
at different doping levels with a special emphasis on the 
possibility of phase separation. Then, we perform the 
similar study for A < 2Jh- After that, we construct the 
phase diagram of the system at A/ J//— doping plane. 

Actually, A, or rather A/t, regularly depends on the 
rare earth radius ta in the series of RC0O3 perovskites 
and increases with decreasing ta- 

III. CHARGE TRANSFER EFFECTS AND 
SPIN-STATE TRANSITIONS: THE CASE OF 
LS-LS GROUND STATE FOR ISOLATED IONS 

Let us first discuss the situation corresponding to 
doped perovskite cobaltites, where Co'^"'" and Co^+ ions 
occupy the sites of a simple cubic lattice. The rela- 
tive number of Co^"*" and Co'^"'" is respectively x and 
1 — X. Let us assume that in the absence of charge 
transfer the cobalt ions of both types are in the LS state 
(A > 3Jh)', this is a typical situation e.g. in the hole- 
doped Lai_a;Sra;Co03, and even more so for smaller rare 
earths R in doped RCoOg^. By promoting some Co^+ 
ions to the IS state, we can have a gain in kinetic energy 
related to the charge transfer from IS Co^+ to LS Co'*"'". 

To treat this situation in more detail, let us introduce 
creation operators aj, and cj, for an electron at the eg 
level and a hole at the t2g level, respectively, at site n 
according to the following rules (choosing Co^^ as the 
vacuum state) 

|0) - |Coi+), e'^---^ ^ Eo, 
al\0) ^ \Co^+), ^ u', 

4|0) = |Coi+), E^^^^E,. (1) 



In terms of these operators, the intermediate-spin state 
of Co'^"'" ions can be constructed in the following way 

\Go^is)^ciai\0), E%+^ =E, + A-Jh = E2. (2) 

Summing up all possible low-energy configurations, we 
can write the following single-site Hamiltonian 

= Eo + {El - Eo)ni + {U' - E^)nl + 

[{E2 - Eo) - {E, - Eo) - {U' - Eo)] <<,(3) 

where n„ = aj,an and = cJ,Cn are the operators de- 
scribing the numbers of electrons at in eg levels and holes 
at t2g levels, respectively. Writing ([3]) in a more compact 
form, we have 

= [Eo + iE,~Eo){n'^-n''J] + 

(A - JH)nl + ;7<(1 - 0, (4) 

where U — U' + Ei — A + Jh ■ Taking the sum over all 
lattice sites and introducing the intersite hopping terms, 
we get 

H = Y,[Eo + {El -Eo- y){ni - n^)] + 

n 

+Ai^< + C/5]<(l-<) 

n n 

-t ^ (aJjOm + h.c) , (5) 

(nm) 

where Ai = A — Jh- 

In Hamiltonian (O, we took into account only the most 
significant hopping integral t describing the transitions 
of electrons from the occupied eg level of IS Co^+ to the 
empty eg level of LS Co^"*". Moreover, we assume that an 
electron can move only without changing the z projection 
of its spin, and, because of the Hund's rule coupling the 
spin of an itinerant {eg) electron and the total spin of 
core {t2g) electrons are parallel to each other. Hence, we 
can assume a ferromagnetic ground state and omit a spin 
index of electron operators. We also neglect the possible 
complications related to the orbital degeneracy of the eg 
level occupied by a single electron; these are not crucial 
for the present problem. 

Such simplified Hamiltonian ^ is quite similar to that 
of the Falicov-Kimball model^^. In model ([5]), we have, 
in fact, an interplay between the electron localization in 
the LS state and the itinerancy in the IS state. This kind 
of interplay was analyzed in detail both analytically^ and 
numerically^^, and a tendency for a nanoscale phase sepa- 
ration was demonstrated. The local (atomic scale) charge 
and spin inhomogeneities related to electronic phase sep- 
aration were also found recently in exact calculations for 
small clusters^. Here, following the technique suggested 
in Refs. ITIfsol we address the specific features of the sys- 
tems with the spin-state transitions. 

The average numbers of eg electrons and t2g holes per 
site (n^) = n'^ and {rJ^) = nP" obey the evident relation- 
ship n'^ —n'^ = X (by electrons we mean here not the real 
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extra electrons, which would create the state Co^"*", but 
the electrons in the initially empty Cg levels, promoted 
there by the LS-IS transition, i.e. we still are dealing 
with the "mixture" of 00"^+ and Co'*+ in the hole-doped 
system). 

Then, the energy per site can be written as 

=£;o(l-x) + £;ix + (i?i}/7V,- (6) 

where 

n n 

-t {a)^am + h.c) . (7) 

(nm) 

Let us consider a homogeneous state corresponding 
to a certain density n*^ of electrons promoted to the 
IS Co'^^ state. We calculate the energy spectrum us- 
ing the Hubbard I decoupling^ in equation of motion 
for the one-electron Green function G"^(n, Hq; t — to) = 
—i{Tan{t)al^^^{to)) for these promoted e„ electrons (anal- 
ogous to the band (b) electrons in Refs. |ll andlssi). In the 
frequency- momentum representation 

^e., ^ ^ u; + ^l ~ Ai - Un'' 

^ (c. + M-Ai-i^i(k))(^ + M-Ai-S2(k)) ' 

(8) 

where 

(9) 

and £(k) is the energy spectrum at ?7 = 0. We choose 
e(k) in the simplest tight-binding form, ignoring possible 
orbital effects and not taking into account the specific 
features of the hopping integrals of eg electrons, [e(k) = 
—2t{coskx + cosfcj, -I- cosfc^) for simple cubic lattice]. 

Using Eq. ^ for the Green function, we calculate the 
densities n'^ and n'^ of Cg electrons and holes. Based on 
these results, we can determine the dependence of the 
total energy on the doping level x. These calculations 
are similar to those performed in Ref. Ill The b electrons 
in Ref. [l| correspond to our eg electrons at IS Co''+ ions, 
whereas the number of localized I electrons in Ref. [l|, 
ni, corresponds to the number of LS Co'^"'" ions, 1 — n^. 
Using this similarity, we could make a direct mapping 
between the two systems. However, in the systems with 
spin-state transitions there exists another homogeneous 
state in addition to that considered in Ref. 1., Namely, 
it corresponds to all Co'^+ ions in an intermediate spin 
state (n'* — 1, n'^ — 1 — x). Formally, in terms of a 
conduction band and localized level, this state can be 
treated as a combination of empty localized level lying 
below the Fermi level of the partially filled conduction 
band, which is in general not possible. In our case, the 
state corresponding to localized level disappears in the 
absence of LS Co^+ ions. 
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Figure 1: (Color online) The energies of type 1 and 2 states 
as function of doping x. The red solid curve corresponds to 
the more favorable in energy homogeneous state, whereas the 
states with higher energies are shown by blue dashed lines 
(see the text). The green dot-dashed line corresponds to the 
energy of an inhomogeneous state, obtained by Maxwell con- 
struction. Ai/zt = 0.2. 



Let us denote the state similar to that in Ref. [J (that 
is the state with coexisting LS and IS Co'^"'" ions) as a 
type 1 state, and the state without LS Co'^"'" as type 2 
state. The energies of these two states as function of 
doping X are shown in Fig.[T]at Ai/zt = 0.2 (z = 6 is the 
number of nearest neighbors) . The type 2 state becomes 
favorable at x > X2- Note that at x > x^ both states, 1 
and 2, are equivalent. We can see that at x < xi there 
are no electrons promoted to the eg level {n^ — 0), see 

Fig.m 

At x > xi the number of eg electrons gradually grows. 
In the absence of type 2 state, this growth would con- 
tinue up to X = X3 when all Co^+ ions would turn to the 
intermediate spin state. At x = a;2, however, the type 
2 state becomes favorable in energy, and the jump-like 
transition to this state occurs. The ground state energy 
E for the homogeneous system as function of x is shown 
in Fig. [T] by red solid curve. At the same time, it is clear 
from this figure that in the doping range < a; < xs 
the inhomogeneous state, being a mixture of states with 
n'' = 1 and n'^ = 0, is more favorable. The energy of this 
mixed state is shown in Fig. [T] by the green dot-dashed 
line. 

The densities of 00"^+ ions in intermediate-spin (n'^) 
and low-spin {1 — n^) states as a function of x are shown 
in Fig. [2] by blue dot-dashed and red solid curves, respec- 
tively. We see a jump- like increase of n*^ at X2, when the 
homogeneous type 2 state would become favorable. The 
behavior of and 1 — n'* at x > a;2 in the absence of 
type 2 state are shown by thin blue and red dashed lines, 
respectively. 

At a small band filling, n'^ ^ 1, we can write an ap- 
proximate explicit expression for the total energy E as- 
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Figure 2: (Color online) The densities of Co ions in 
intermediate- {rf — njg co3+) s^nd low-spin [1 — — 
riLg co3+) homogeneous states as a function of x aX A\/zt = 
0.2. At xi < X < X2, both spin states of Co"^""" coexist, whereas 
at X — X2, there occurs a jump-like transition to the purely 
IS state of Co^+ ions. Dashed lines illustrate the possible be- 
havior of Tijg co3+ ^i^d rija^ co3+ for the type I state similar 
to that described in Ref. [J 



values. Note here that the values of magnetic moment 
under discussion correspond to macroscopic phase sepa- 
ration, that is the characteristic sizes of inhomogeneities 
are much larger than the lattice constant. It is indeed so 
at relatively large x (exceeding the percolation threshold 
for the phase with itinerant charge carriers). At small x, 
it is naturally to expect that the phase-separated system 
will consist of small droplets (spin-state polarons) con- 
taining only one Co^"*" ion surrounded by IS Co'^"*'. In 
the latter case, the magnetic moment per Co^"*" should 
be larger than that corresponding to the macroscopic 
phase separation. This could be the case for spin po- 
larons in low-doped Lai_a;Sra;Co03 observed in Ref. 1271. 
where the polarons with the magnetic moment equal to 
13/is seem to be the most probable. From our consider- 
ations, one should expect that the value of the moment 
per of Co*+ should become smaller with the increase of 
doping X. The exact calculations for small clusters also 
demonstrate that in a suitable range of parameters the 
saturated magnetic moment can exist at relatively low 
temperatures also in atomic-size doped clusters of vari- 
ous geometries^. 



suming that the Fermi surface is spherical 

E-Am^^ tzn^n^ + - (367r V) (n^f'' . (10) 
5 

The density n'^ of itinerant electrons is determined by 
minimization of Eq. (jlOp with respect to n'^ taking into 
account that = x + n'^. It can be easily shown that 
the solution for the energy minimum corresponding to 
rf ^ Q can exist only if Ai < tzx. This means that at 
Ai/iz > 1 the LS Co^+ ions can not be promoted to the 
IS state at any doping x. 

The dependence of n'^ on doping x determines the be- 
havior of magnetic moment of Co ions. Indeed, the LS 
Co'^^ ions correspond to zero magnetic moment, S — Q, 
while the doping leads to creation of LS Co^+ ions 
(5 = 1/2) and also provides the promotion of some 00"^+ 
ions to the IS state [S — 1). So, the data presented in 
Fig. [2] could be redrawn in terms of magnetic moment per 
dopant (or, in other words, per 00"*+), see Fig.|31 For the 
homogeneous state, we see in Fig. |3] that the jump-like 
transition in the density n'^ of itinerant electrons mani- 
fests itself in a jump of magnetic moment. At the same 
time, in the phase-separated state, the magnetic moment 
per Co^+ ion remains constant since both the content of 
the phase with IS Co^+ and the number of Co^+ ions are 
proportional to x. The value of magnetic moment per 
Co**^ is determined by the value of x, where the green 
dot-dashed line in Fig. [1] touches the curve corresponding 
to the energy of the homogeneous state. Both the height 
of the jump for the magnetic moment in the homoge- 
neous state and the value of magnetic moment per Co^+ 
in the phase-separated state depend drastically in the pa- 
rameters of the model, especially on the hopping integral 
t. In Fig. [21 we see that the increase in t by a factor 
of two leads to a pronounced growth of both mentioned 
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Figure 3: (Color online) Magnetic moment per one Co*^ ion 
versus doping x at different values of the hopping integral t. 
The red curves correspond to the homogeneous states. The 
behavior M in the phase-separated state is shown by the blue 
dashed line. 
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Thus, we demonstrated that the spin-state transitions 
in hole-doped cobaltites can be described based on the 
model involving the coexistence and competition of lo- 
calized and itinerant electron states. In contrast to the 
similar model for manganitesii^S,, this model allows the 
possibility of a jump-like transition to the purely itiner- 
ant state corresponding in the case of cobaltites to the LS 
— >IS transition for all Co'^+ ions. However, at lower dop- 
ing, before reaching this homogeneous metallic state with 
all Co ions magnetic, the phase-separated state comes 
into play, in which only a part of Co'^"'" ions is promoted 
to the IS state, doped holes being located in these regions. 
Experimental data on Lai_a;Sra;Co032^i^i2ii2^i2^iS2 seem 
to be in agreement with this picture. 

IV. CHARGE TRANSFER EFFECTS AND 
SPIN-STATE TRANSITIONS: THE CASE OF 
HS-HS GROUND STATE FOR ISOLATED IONS 

Let us now discuss the situation at A < 2 J//, when in 
the absence of electron hopping it is favorable for both 
Co'^"'" and Co^+ to be in the HS state. The charge transfer 
becomes possible only if we promote a hole to the Cg level 
of Co^+ and transform such an ion from HS to IS state. 

So, in this case, instead of electron hopping from IS 
Co'^"'" to LS Co^+ , we have the electron hopping from the 
HS Co3+ to IS Co4+, or the hole hopping from IS Co^+ to 
HS Co'^+ (this representation is more convenient here). 
Using this analogy, we can choose the HS state of Co^+ 
as a new vacuum state and write relationships similar to 
(H]) and © as 

|0) = |Cot,+s), -B^™^) - -f 2A - 6Jh = ^0, 

gt^io) = |Co5+), £;(5+) = ii', 

dlcl\0) = \Coj+),E^/+^ ^Eo + A-2Jh^E2. (11) 

The corresponding single-site Hamiltonian can be 
found by the following substitution in ([3]) and ([1]): 

Eq, El, E2,U Eq, El, E2,U and also "n^^-n ~* 

As a result, we can rewrite the Hamiltonian ([5]) in the 
following form 

H = J2[Eo + {Eo~El-^l){fi'^^n':,)] + 

n 

n n 

-t J2 (45m + h.c.) . (12) 

(nm) 

Here, A2 = AJh — A is the energy difference between 
the IS and HS Co'*"'' ions, cj,, Cn are creation and annihi- 
lation operators for a hole promoted to the Cg level of IS 
Co"*"*" at site n, = cJ,Cn, and = aj,a„ is the opera- 
tor describing the number (0 or 1) of additional localized 



t2g electrons at site n (aj,, an are creation and annihila- 
tion operators for such electrons). The average numbers 
of electrons and holes per site obey now the relationship 
fi^ — h'^ = 1 — X. 

In this case, the energy per site (|6]) can be rewritten 

as 

E^^^^ ^Eo{l-x)+Eix+{H2)/N, (13) 

where 

H2 = J2 (2^ - ^■^H + 2Jh{K - "n)) + 

n 

II n 

-t (4am + h.c.) . (14) 

(nm) 

Note that the difference i?^^) — E^^\ does not depend 
on the choice of Eq and Ei, this fact will be helpful in 
constructing the phase diagrams in the next section. 

Thus, the behavior of the system energy and charge 
carrier densities, h'^ and n'^ are similar to those shown in 
Figs.[T]and[21 In these figures we should replace n'^ — > n'^, 
nf' ^ Tip- , and a; — » 1 — x, that is, the densities of Co"^"*" 
ions in IS (n"^ = njg ^0^+) ^-i^d LS (1 — n'' = Co-''+) 
states become here the densities of Co'*+ ions in IS [n^ = 
"-IS, Ga-^+) HS {1 — ff = riHs, Co*+) states. Note also 
that such an exact similarity between the LS-LS and HS- 
HS cases appears since we, in fact, deal with the spinless 
fermions (the spins of charge carriers are parallel). So, 
we have an electron-hole symmetry between an empty eg 
level at LS Co'^"'" and a completely occupied such level at 
HS Co3+. 



V. PHASE DIAGRAMS 

Based on the results of the previous sections, we can 
summarize the behavior of the system as function of dop- 
ing at different values of the A / Jh ratio and to draw the 
corresponding phase diagram. The form of this phase di- 
agram depends drastically on the characteristic values of 
the hopping integral t. The general features of the evo- 
lution of the system with doping from one homogeneous 
state to another are illustrated in Fig.[31 At rather small 
t {t/ Jh ^ 1), see Fig. |4^, we have clearly defined re- 
gions of the phase diagram corresponding to A > 3 J// 
and A < 2Jh (corresponding to the situations discussed 
in sections IlIII and HVl respectively). In each of these re- 
gions, the variation of doping leads to the transitions be- 
tween the phase with only localized carriers to the phase 
when some charge carriers are delocalized and, eventu- 
ally, to the phase when all charge carriers are itinerant. 
These two regions, with A > 3 J// and A < 2 J//, are sep- 
arated by the phase with Co'^"'" in LS {S — 0) and Co"*"*" 
in HS (5* = 5/2) states, with the charge carriers localized 
because of the spin blockade^*. At larger t (tj Jh ^ 1), 
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see Fig. |4]d, the latter intermediate region collapses at a 
certain doping range, and a direct spin-state transition 
between the phases with fully delocalized charge carriers 
becomes possible. 



< 





/ / 


■ no itinerant charge carriers (ChC) 






some ChC arc itincfant 




/ 

/ 


all ChC are itinerant 


/ 


",,sc„'=0 


localized ChC, COj^^'/Co^^"' 


all ChC are itmerant 


■ / 

/ 




/ 

/ 


some ChC are itiner^ilt 

,/^'co,;-/co,;Vcp;;^ 


no itinerant ChC 






/ 1 . / 1 





0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



(a) tlJi-_ 



no itinerant charge carriers (jChC) 
Co, '7Co, 



^omc ChC gre itinerant 



all ChC are itinerant 




all ChC are itinerant 
n .-0 



some ChC'are itinerant' 



rs HS' 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



(b) t/JH = 1.5 



Figure 4: (Color online) Possible homogeneous states of the 
system under study at different values of the hopping integral 
t. The boundaries between analogous phases in the upper and 
lower parts of the phase diagram are shown by the same lines 
(solid, dashed, or dot-and-dash). 

The form of the phase diagram changes if we take into 
account the possibility of phase separation. The cor- 
responding phase diagrams drawn at different values of 
tj Jh are shown in Fig. [51 We see that instead of phases 
with partially and fully delocalized charge carriers, there 
appears a broad regions of phase separation where the do- 
mains of fully localized and fully delocalized charge car- 
riers are intermixed. Again, at rather small t {t/ Jh ^1), 
we have and intermediate region where the charge carries 
are localized at any doping level (with Co'^"'" and Co^+ 
in LS and HS states, respectively). This intermediate re- 
gion gradually disappears with the growth of the hopping 



< 
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PS 1 



localized charge carriers (ChC), Co^^ ^^'^hs 



PS n 



0.0 



0.2 



0.4 
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(a) t/J„ = l 



PS I 




PS II 



(b) t/Ji, 



0.6 



1.3 




(c) t/Jn = 1.5 

Figure 5; (Color online) Phase diagrams including the phase- 
separated states of the system under study at different values 
of hopping integral t. PS I is the phase-separated state includ- 
ing the regions without itinerant charge carriers, correspond- 
ing to LS Co^"*" , and those with completely delocalized charge 
carries promoted to IS Co^"*". PS II is the similar phase- 
separated state where the regions with and without itinerant 
charge carriers correspond to IS Co'*"'' and HS Co*+, respec- 
tively. Regions 1 at panel (b) correspond to the charge carriers 
located at HS Co*^ and LS Co^^ where the charge transfer 
between Co sites is suppressed due to the spin blockade. 
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integral t. 

Let us note here that the phase diagram along the A 
axis could be reproduced varying the average ionic radius 
of the rare-earth ions in cobaltites (see, e.g. Refs.lSTllssh. 

Note also that the long-range Coulomb interaction re- 
lated to the charge disproportionalization in the phase- 
separated state can reduce the doping range of the phase 
separation and modify the form of the phase diagram 
shown in Fig. [5l 



dependent on the ratio of of the electron hopping integral 
t and and the Hund's rule coupling constant Jh- 

Here, we did not analyzed in detail the possible struc- 
ture of the phase-separated state. However, for the corre- 
sponding model describing doped manganites, the calcu- 
lations^ and numerical simulations^ taking into account 
the surface and long-range Coulomb contributions to the 
total energy lead to the characteristic size of nanoscale 
inhomogeneities of the order of several lattice constants. 



VI. CONCLUSIONS 
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